Wave number dependence of the transitions between travelling 
and standing vortex waves and their mixed states in the 

Taylor-Couette system 

A. Pinter*, M. Liicke, and Ch. Hoffmann 

Institut fiir Theoretische Physik, Universitdt des Saarlandes, 
Postfach 151150, D-66O4I Saarbriicken, Germany 
(Dated: Jufy 18, 2008) 

Abstract 

Previous numerical investigations of the stability and bifurcation properties of different nonlin- 
ear combination structures of spiral vortices in a counter-rotating Taylor-Couette system that were 
done for fixed axial wavelengths are supplemented by exploring the dependence of the vortex phe- 
nomena waves on their wavelength. This yields information about the experimental and numerical 
accessability of the various bifurcation scenarios. Also backwards bifurcating standing waves with 
oscillating amplitudes of the constituent traveling waves are found. 
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Recently the stability exchange between travelling waves (TWs) |l| and standing waves 
(SWs) of spiral vortices in the Taylor-Couette system has been investigated by full numer- 
ical simulations and a coupled amplitude equation approximation j^. TWs and SWs have 
a common onset as a result of a primary, symmetry degenerate oscillatory bifurcation. The 
SW solution is a nonlinear superposition of mirror symmetric, oppositely propagating TWs 
with equal amplitudes. At onset either the TW or the SW solution is stable 3|, Then, 
at larger driving there is a secondary bifurcation that leads to a stability exchange between 
the two solutions. This exchange is mediated by mixed patterns that establish in solution 
space a connection between a pure TW and a pure SW. The mixed structures consist of 
a superposition of oppositely propagating TWs with temporally constant, nan- equal ampli- 
tudes. The TWs investigated in [2] are initially stable while the SWs gain stability later 
on. 

There is a second variety of mixed states in which the TW amplitudes oscillate in time in 

counterphase. This stable solution bifurcates out of the SW at even higher driving rates via 

n 

a Hopf bifurcation [5] in which the aforementioned SWs lose their stability. These results 
lave been found by full numerical simulations of the vortex flow in a Taylor-Couette system 

1 Q 

6, |7| with counter-rotating cylinders of radius ratio rj = 0.5 with methods described in 
8|. The calculations were done for a fixed axial wavelength A by imposing axially periodic 
boundary conditions 9|. 

Here we investigate and show how stability, bifurcation properties, and the spatiotem- 
poral behavior of the aforementioned structures change with A. Thus, this report provides 
information for future simulations and experiments with finite length set-ups and, say, non- 
rotating lids that close the annular gap between the cylinders at the ends: Since the height 
of the system influences the wavelength of the vortex structures and with it their properties 
the prior knowledge of their A-dependence is of significant interest. 

Structures — The following structures have been investigated: (i) Forward bifurcating 
TWs consisting of left handed spiral vortices (L-SPI) or of right handed spiral vortices (R- 
SPI) that are mirror images of each other. L-SPI (R-SPI) travel in the annulus between the 
two cylinders axially into (opposite to) the direction of the rotation frequency vector of the 
inner one, i.e., in our notation upwards (downwards) [8]. (ii) Forward bifurcating SWs that 



consist of an equal-amplitude nonlinear combination o 
called ribbons (RIBs) in the Taylor-Couette literature 



L-SPI and R-SPI. These SWs are 



lo|,[iil. 



[Hi) So-called cross-spirals 
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(CR-SPI), i.e., combinations of L-SPI and R-SPI with different stationary amplitudes. They 
provide a stabihty transferring connection between TW and SW solution branches [2], [gj. 
And, finally, {iv) oscillating cross-spirals (0-CR-SPI). Therein, the amplitudes of the TW 
constituents of the SW, i.e., the amplitudes of L-SPI and the R-SPI oscillate in counterphase 
around a common mean [5|. The vortex structures (z)-(ro) are axially and azimuthally 
periodic with axial wave number k = 2tt/X and azimuthal wave number M = 2 in our case. 

Control- and order parameters — The control parameters are the Reynolds numbers 
i?i > and R2 < defined by the rotational velocities of the inner and outer cylinder, 
respectively. As order parameters we use the amplitudes 

A{t) = U2,iit), Bit) = n2,-i(t) (1) 

of the dominant critical modes of the radial velocity u at mid-gap in the double Fourier 
decomposition in azimuthal and axial direction. In Eq. ([T]) the indices m = 2 and n = ±1 
identify azimuthal and axial modes, respectively. Note that for SPI, CR-SPI, and RIB 
structures investigated here the moduli in Eq. are constant. On the other hand, in 
0-CR-SPI the moduli \A{t)\ and oscillate in counterphase around a common mean. 

Therefor we use the difference of the squared moduli D{t) = (|y4(t)p — \B{t)\'^)/2 and its 
oscillation amplitude D to describe the bifurcation from the RIB solution {D = 0) to O- 
CR-SPI (5^0). 

X-dependence of the bifurcation scenario — Fig. [1] shows the A-dependence of the bifur- 
cation thresholds Ri for M = 1 and M = 2 SPI and RIB and M = Taylor vortex flow. 



12 1 of the basic circular Couette 



These results were obtained from a linear stability analysis 
flow. For the two characteristic values R2 = —540 and —605 shown there the M = 2 SPI 
and RIB have the lowest threshold for a wide range of 0.8 < A < 2.1. 

In Fig. [2] we show for different A bifurcation diagrams as functions of R2 for a fixed 
Ri = 240. For 1 < A < 1.2, figure parts (a)-(c), the bifurcation properties are quite 
similar: SPI and RIB bifurcate supercritically out of the circular Couette flow, SPI (RIB) 
are unstable (stable) at onset, and there is no stability exchange in the range of R2 of Fig. 

12 

By contrast, for A = 1.3 and 1.4, in (d) and (e), respectively, there are different interesting 
stability exchanges. Here, L-SPI (A ^ 0, B = 0) and R-SPI {A = 0, B ^ 0) are initially 
stable while the RIB state {A = B) is initially unstable. But RIB gain stability almost 
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immediately thereafter: the stabihty transfer from L-SPI or R-SPI to RIB is mediated within 
a very small interval by the L-CR-SPl {\A\ > \B\) or the R-CR-SPl {\B\ > \A\) solution, 
respectively [2]. For larger driving, the RIB lose stability again, when stable oscillating 
structures, 0-CR-SPI, appear via a Hopf bifurcation Jsl. Note, however, that 0-CR-SPI 
bifurcate forwards for A = 1.3 but backwards for A = 1.4, cf. details further below. 

In Fig. [3] we show bifurcation diagrams as a function of A for Ri = 240 and two different 
i?2 indicated by arrows in Fig. [2l In the case of R2 = —595 (right arrow in Fig. [2]) SPI 
are unstable and RIB are stable for all A. For R2 = —605 (left arrow in Fig. [2]), on the 
other hand, this stability situation — RIB are stable and SPI are unstable — applies only 
as long as A < 1.25: Then, with increasing A, a stability exchange between RIB and SPI via 
CR-SPI occurs that is reflected also at the very left end of the bifurcation diagram in Fig. [2] 
(d). So, the interesting stability exchange between TWs and SWs occurs in a rather narrow 
wave number band around A ~ 1.3. 

In figures |2] and [HI we showed the case where stability is transferred from L-SPI to RIB. 
The symmetry degenerated situation where stability is transferred from R-SPI to RIB via 
R-CR-SPI is obtained by exchanging (the symbols for) A and B in these figures. 

Phase diagram for A = 1.3 — In view of the above discussed stability exchange process 
we take a more detailed look at the case A = 1.3 for which previous calculations have been 

m 



done only at the two Reynolds numbers -Ri = 200 and 240 



13j . To that end we provide 



in Fig. m the phase diagram of the stable, aforementioned M = 2 vortex structures with 
fixed A = 1.3 in the R1 — R2— parameter plane. Stable M = 2 SPI appear first via a primary 
forwards bifurcation at the lower left border of the red stripe in Fig. H] [l^ . Then, for a fixed 
Ri ^ 190, we have observed with increasing R2 always the same stability transfer sequence: 

SPI ^ CR-SPI ^ RIB ^ 0-CR-SPI. 

For lower Ri, however, the existence range of stable M = 2 structures seems to be 
more and more confined from above by the appearance of M = 1 modes at the respective 
bifurcation threshold (dashed line): With decreasing Ri first the 0-CR-SPI and then the 
RIB and CR-SPI areas are pinched off successively. 

Note that in all cases the CR-SPI stripe is extremely thin (cf. the blow-up bar in Fig. 
S]) whereas the 0-CR-SPI area being quite large should facilitate a respective experimental 
observation. 
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Backwards bifurcating 0-CR-SPI — As noted already in the discussion related to 
Fig. [2](e), we have found for the first time backwards bifurcating 0-CR-SPI. In Fig. [5] we 
display for fixed A = 1.4 and R2 = —605 as a representative example the bifurcation prop- 
erties of this new scenario. Fig. [5](a) shows the squared moduli |Ap and as a function 
of -Ri. Diamonds show stable RIB. They have obtained their stability from the SPI via a 
CR-SPI branch connection at smaller Ri outside the plot range of Fig. [51 

The + and — signs denote the maximal and minimal amplitudes, respectively, of the 
modes A and B that oscillate in counterphase in the 0-CR-SPI. The hysteresis in the 
transition between stable RIB and stable 0-CR-SPI is best visible with the order parameter 
D in Fig.EKb). 

Conclusion — Our results show that the mixed states of stationary CR-SPI and of 0-CR- 
SPI should be observable in experimental setups or in finite lengths numerical simulations 
when the wavelength of these vortex structures lies in the interval of 1-3 < A < 1.4. Therein 
0-CR-SPI are stable in a wide range of control parameters. They bifurcate either forwards 
or, as we have found here, backwards out of the RIB state of standing waves. CR-SPI 
solutions, on the other hand, exist only in a rather small interval of control parameters. 
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FIG. 1: Bifurcation thresholds i?^ of M = Taylor vortex flow and of SPI and RIB with azimuthal 
wave numbers M = 1 and M = 2 versus axial wavelength A for different R2 as indicated. 
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FIG. 2: (Color online). Influence of the axial wavelength A on the bifurcation behavior for fixed 
Ri = 240. Solid (open) symbols denote stable (unstable) structures. Symbols show L-SPI (red 
circles, vl / 0,5 = 0), L-CR-SPI (violet triangles, \A\ > \B\), RIB (blue diamonds, A = B), and 
stable 0-CR-SPI (magenta). In the latter |^(t)| and |-B(t)| oscillate between the + and — symbols. 
In (d) and (e) only stable structures are shown. The arrows refer to the R2 values of Fig. [3l 
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FIG. 3: (Color online). Bifurcation diagrams of L-SPI (red circles, A 7^ 0, i? = 0), RIB (blue 
diamonds, A = B), and L-CR-SPI (violet triangles, \A\ > \B\) as a function of the axial wavelength 
A for fixed Ri = 240 and the two values of R2 that are indicated by arrows in Fig. [2j Solid (open) 
symbols denote calculated stable (unstable) structures. Full (dashed) lines for stable (unstable) 
solutions branches result from a spline interpolation. However, the violet L-CR-SPI branches are 
schematic. 
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FIG. 4: (Color). Phasediagram of stable SPI, CR-SPI, RIB, and 0-CR-SPI. In the colored areas 
these vortex structures (as indicated by the color code) with azimuthal wave number M = 2 and 
axial wavelength A = 1.3 are stable. The white area has not been investigated in this work. The 
basic state of circular Couette flow is stable in the dotted region: Lines are marginal stability 
boundaries of circular Couette flow against growth of vortex flow with A = 1.3 and azimuthal wave 
numbers M = 2 (full), M = 1 (dashed), and M = (dotted). The bar shows a blow-up of the 
region between the black arrows. 
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FIG. 5: (Color online). Backwards bifurcation of 0-CR-SPI (magenta circles and + and — signs) 
out of RIB (blue squares) for A = 1.4 and R2 = —605 as a function of Ri. In (a) squared 
amplitudes \A\'^, \B\'^ are shown with symbols as explained in Fig. [2l The oscillation amplitude D 
oi D = (|Ap - |Bp) /2 is shown in (b). For RIB D = D = 0. Full (dashed) lines denoting stable 
(unstable) solutions were obtained by spline interpolation. 



